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^ ; Abstract 

o 

OO ' Axions are produced during a period of dilaton-driven inflation by amplifi- 

cation of quantum fluctuations. We show that for some range of string cos- 
mology parameters and some range of axion masses, primordial axions may 
^ , constitute a large fraction of the present energy density in the universe in the 

Q I form of cold dark matter. Due to the periodic nature of the axion potential 

' energy density fluctuations are strongly suppressed. The spectrum of primor- 

dial axions is not thermal, allowing a small fraction of the axions to remain 
^ . relativistic until quite late. 



1 



Typeset using REVTgX 



I. INTRODUCTION 



Axions are hypothetical particles, invented to solve the strong CP problem Many 
ongoing experimental efforts aim to detect axions, but have yet to produce evidence for their 
existence. Cosmological and astrophysical implications of axions are well studied 0,^, in 
particular, axions are among the leading candidates for providing the missing dark mass in 
the universe. 

String theory possesses many axion candidates of the "invisible axion" type [0]. We 
will focus on the "model-independent axion" ||^ which respects a Peccei-Quinn symmetry to 
all orders in perturbation theory. There are no good general arguments that QCD provides 
the dominant contribution to the potential energy of any of the stringy axions, including 
the model independent axion. However, there are some theoretical conditions under which 
the model-independent axion could be the axion that solves the strong CP problem. In any 
case, we will assume that this is indeed so. 

We consider axion production in models of string cosmology which realize the pre-big- 
bang scenario |^,^. In this scenario the evolution of the universe starts from a state of 
very small curvature and coupling and then undergoes a long phase of dilaton-driven kinetic 
inflation and at some later time joins smoothly standard radiation dominated cosmological 
evolution, thus giving rise to a singularity free inflationary cosmology. Axions are produced 
during the period of dilaton-driven inflation by the standard mechanism of amplification of 
quantum fluctuations [0. 

The spectrum of relic axions depends on their potential and interactions, and on some of 
the string cosmology model parameters. By applying simple constraints, such as requiring 
that the energy density of the universe does not exceed the critical density at different 
stages of the evolution, we are able to constrain parameters of string cosmology models and 
axion potential and find a consistent parameter range in which most of the energy of the 
universe today is in the form of cold dark matter axions whose origin is quantum fluctuations 
from the pre-big-bang. This consistent range overlaps with the range in which relic gravity 
wave background produced during the dilaton-driven inflationary phase could be detected 



by planned gravity wave experiments |[TT|]. The same parameter range could perhaps lead 



to formation of observable primordial black holes [O . 



II. THE MODEL 

We assume that the model independent axion receives the dominant part of its potential 
from QCD instantons, therefore, roughly speaking, the axion is massless until the universe 
cools down to a temperature Tqcd ~ ^qcd at time t = tqcn- A more sophisticated 
estimate as in will be used later. The axion then develops a periodic potential of 
overall approximate strength Ag^^ (recall that Aqcd ~ 200MeV), and period which is 
apriori a free parameter /pg. The scale fpq is typically less than lO^^GeV resulting in 
an axion mass rua ^ 10~^^eV. Astrophysical constraints further bound fpq from below 
fpQ ~ lO^GeV resulting in an axion mass rria ^ lO'^eV. We will discuss cosmological 
constraints on fpq in more detail later. The explicit form of the axion potential that we will 
assume V{tlj) = |Vo (l — cos(^)^, depends on two parameters, ipo, related to the Peccei- 
Quinn scale fpq, and Vq = mlipQ. For QCD axions Vq ~ f^^l ~ ^qcD- 
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The model of background evolution we adopt in this paper is a simplified model used 
also in The evolution of the universe is divided into four distinct phases, the first phase 



is a long the dilaton-driven infiationary phase, the second phase is a high- curvature string 
phase of otherwise unknown properties, followed by ordinary Friedman-Robertson- Walker 
(FRW) radiation dominated (RD) evolution and then a standard FRW matter dominated 
(MD) evolution. We assume throughout an isotropic and homogeneous four dimensional 
flat universe, described by a FRW metric. The model is described in full detail in |TB|, and 



we reproduce here only its important features. We note in particular that the axion field is 
assumed to have a trivial vacuum expectation value during the infiationary phase. 

The dilaton-driven infiationary phase lasts until time t = tg, and while it lasts the 
scale factor a{t) and the dilaton are given by the solution of the lowest order string- 
dilaton-gravity equations of motion, the so-called branch vacuum. The string coupling 
parameter e'^ = g stringy the models that we consider is, of course, time dependent. Both 
curvature and coupling e'^ are growing in this phase, which is expected to last until curvatures 
reach the string scale and the background solution starts to deviate substantially from the 



lowest order solution. For ideas about how this may come about see [0 . 

The string phase lasts while tg < t < ti. We assume that curvature stays high during 
the string phase. As in |[T^, we assume that the string phase ends when curvature reaches 
the string scale Mg, H{ti) ~ M^. We parametrized our ignorance about the string phase 
background, as in by the ratios of the scale factor and the string coupling g{t) = e'^'^*^/^, 
at the beginning and end of the string phase zs = ai/as and gi/gs, where gi = 6"^^*^^/^ 
and gs = e'^^*'^^/^, where as = a{ts) and (ps = 0(^s)- We take the parameters to be in a 
range we consider reasonable. For example, zs could be in the range 1 < < e^^ ~ 10^°, 
to allow a large part of the observed universe to originate in the dilaton-driven phase, and 
gi/gs > 1; assuming that the coupling continues to increase during the string phase and 
10^'^ < (71 < 10^^ to agree with the expected range of string mass (see e.g. ||15[|). Some other 
useful quantities that we will need are ui, the frequency today, corresponding to the end of 
the string phase, estimated in to be ui ~ lO^^Hz, and the frequency us = uji/zs, the 
frequency today corresponding to the end of the dilaton-driven phase. 

Standard RD phase and then MD phase are assumed to follow the string phase. The 
dilaton is taken to be strictly constant, frozen at its value today. 

We have presented our assumptions about background evolution and axion potential in 
great detail, and will use them as presented, even though many of the assumptions can be 
either relaxed (without affecting dramatically our results), or improved to take into account 
additional known effects. However, each change adds an additional level of complication by 
adding parameters and assumptions, and we preferred to keep the discussion as simple as 
possible to capture the essential physics. Nevertheless, we do mention from time to time a 
possible alternative or generalization. 

III. PRIMORDIAL AXIONS 

The spectrum of axionic perturbations produced during the dilaton-driven inflation is 
approximately given by |TB[ (see also |TB|,[T^) 
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f k \ 

ipk = Mkas^gs^ks^'"^ ( — ) ,k < ks, (3.1) 



where Afk = 2^^^T{y/3)/ and ks = ujsa,(t). The r.m.s amphtude of the perturbation 
in a logarithmic k interval is defined in a standard way 6-ipk = k^^'^'4'k- Using the relation 
ks/o-s = ki/ai = uji{ti), the assumption uji{ti) = Mg and the relation between the string 
mass and the Planck mass in weakly coupled string theory, Mg = MpQi we obtain 

5^l^k = MkM, -(^] , (3.2) 

9s \ksj 

for perturbations outside the horizon, and 

*V..=Ar.M,£lfA)""'^(^), (3^3) 

for perturbations that have reentered the horizon during RD, before the axion potential is 
generated. 

The ratio of energy density in axions per logarithmic frequency interval to the critical 
density -j^, for perturbations that reenter the horizon during RD, before the axion potential 
is generated, is given by ||13[ 



'^^1 — — ' w < ws, (3.4) 



dhiu \gs J \uJsJ 

where C is a numerical factor which we will ignore in the following. Note that the spectral 
index 3 - 2^/3 ^ -0.46 is negative, and therefore most of the energy is contained in the 
low-frequency modes. 

The total energy density within the horizon, at a given time, is dominated by the lowest 
frequency which is just reentering the horizon, 

■^sW / X 2 / rr/.N \ 3-2V3 

H{t) 

where us(t) = uji(t)/zs. If H{t) > usit) then the total energy density in axions produced 
during dilaton-driven phase simply vanishes (We will discuss an estimate for the axions 
produced during the string phase later on). Since H{t) oc T^(t) and uJs{t) oc T{t), and since 
H{ti) = Ms, then ^ = T^/{Tzs), and therefore 



2 / , , X 2V3-3 



2i9i\ / Ms 



To ensure standard RD cosmology at late times we must require that the energy den- 
sity in axions remains smaller than critical f2a < 1. It is enough to require this at the 
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lowest temperature possible, i.e., at the temperature just as the axion potential is gener- 
ated, naitgcD) ^ gl { tqcdzs T ^ ^ < 1- Usiug M, = QiMp, we obtain the following 
condition 

9iMp 



zs > 



QCD 




l/(2v^-3) 

(3.7) 



Note that since ~ 10^^, unless zs is large enough condition ( p.7|) can be satisfied 

only if gi is unacceptably small. There is, however, a reasonable range of parameters for 
which condition ( |3.7| ) is indeed satisfied, for example, if gi = 10~^, gs/gi = 1/10, then ( p.7| ) 
implies > 3 x 10^ and if gs — gi — 0.01 then > 3 x 10^. If gi ~ gs, condition 
( p.7[ ) simplifies to zs ^ g\Mp/TQCD- When ( |3.7| ) is saturated, axions provide near closure 
density of the universe just before the axion potential is generated. Condition ( p.7| ) is valid 
for standard adiabatic RD evolution. If some intermediate period of matter domination or 
entropy production is assumed, condition ( p^ ) is relaxed. 

So far we have considered only axions that were produced during the dilaton-driven phase 
and ignored axions that were produced during the subsequent string phase. We would like 
to show that it is reasonable to neglect axion production during the string phase by giving 
an estimate based on the extrapolation used in |jl^ , which assumes constant H and during 
the string phase. The resulting energy density fi^^ is given by 

9i\ — \[ — ] , ujs <uj<uji, (3.8) 



d\n.uo 



where the spectral index C, = ln{gi/gs)/ Itlzs, is positive and therefore the energy density 
decreases with frequency. The total additional energy in axions produced during the string 

phase / ^^dlnuj, is up to a numerical factor ~ gf f ^, which is indeed negligible (for 

large zs and reasonable gi, gs) compared with the energy density (|3.6| ) in axions produced 
during the dilaton-driven phase. The same conclusion is expected as long as the spectrum 
of axions produced during the string phase continues to decrease. 

If the model-independent axion were to remain massless its total energy just before 

matter radiation equality would be given by ila{teq) — gf (^^^ ( t*% ) ' (^^call that 
Teq ~ leV). Axions would overdose the universe and lead to an unacceptable cosmology, 
unless the parameters of string cosmology, and in particular zs are pushed to uncomfortable 
values. 



IV. DARK MATTER AXIONS 

We turn to discuss the effects of the axion potential as it turns on when the universe cools 
down to QCD temperatures. If we try to approximate the axion potential by a quadratic 
potential, a common practice in most investigations, we encounter a puzzle. The axion 
energy density becomes formally divergent as soon as the axion potential turns on! (if we 
assume that the dialton-driven phase lasted only a finite time then the formal divergence is 
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replaced by a singular dependence on the duration of the dilaton-driven phase). The relative 

dlni. 

we reproduce here its low frequency part, 



energy density in axions, assuming a quadratic potential , was computed in |]T3| and 



dlnuj 




LJ<UJra,UJs- (4.1) 



Once the potential is generated, all the low frequencies reenter the horizon at once, so to 

dine 

amplified frequency Umin, which is either zero, if the duration of the dilaton-driven phase is 



obtain the total energy density inside the horizon we need to integrate from the minimal 



infinite, or exponentially small if the duration is finite but large, ^l^^{t) = J "^^"^ dlnu. 



The lower frequency part of the spectrum yields a divergent contribution, proportional to 
t^min^"^^ (recall 3 — 2y/3 ^ —0.46). This result does not make sense. 

The resolution of the puzzle depends crucially on the periodic nature of axion potential 
V{t/j) = |Vo (l — cos{^)^. This point was first understood by Kofman and Linde |jT8|, and 



we have adopted their ideas to our particular situation. First, the total potential energy 
is limited to Vq and does not continue to increase indefinitely as the axion field increases, 
providing a "topological cutoff" on the total axionic energy density and as important, large 
fluctuations in the axion field are also "topologically cutoff' , producing exponentially small 
energy density perturbations. Large fluctuations lead to a uniform distribution of the axion 
field inside the horizon, with very small statistical fluctuations. 

The axion potential is highly non-linear, therefore it is not possible to solve the pertur- 
bation equation mode-by-mode. In [Q, the following strategy is suggested. Consider the 
axion field ijj{f,t) at the time when the axion potential is turned on. The low k Fourier 
modes ipk, k/aitqcD) < H{tQCD), provide an essentially constant field ipc across the horizon. 
The value of ipc is random, and in our case it is determined statistically by a Gaussian 

/ aH 

distribution P{ipc) with zero average and standard deviation = w / dlnklSipkl"^. Since 

in all cases that we will be interested in, ipQ < lO^^Gel^, and 6tpk > Mp for all k < aH (see 
eq.(|3.2|)), the width cTc is much larger than the period of the axion potential Uc ^ ^nipQ. 
The constant value ipc becomes essentially uniformly distributed among all possible values. 
The average energy density in the non-relativistic part of the axion field pa = {V {tp{r))) is 
given by 



Pa = J IVo (^1 - COS (^|-^ ^ P^d^ ~ IVo, (4.2) 



with exponentially small corrections. Note that eq.( [4.2| ) is valid for all reasonable values of 
zs and gs, gi- Regardless of the fraction of relativistic axions which exists at tqcn, the low 
momentum modes with wavelength larger than the horizon contribute a constant energy 
density. The procedure that we outlined above can be repeated for any scale i, separating 
modes of ip{r, t) = ipd^'^) + V'? where ip contains only modes with k < 

The constant axion field ipc starts to coherently oscillate around the nearest minimum 
of the potential. Using completely standard arguments [p|,^,p!9| , we may obtain a bound 
on nia (or equivalently on i/jq) by requiring that the energy density in the coherent axion 
oscillations be subcritical at the beginning of MD epoch. This requirement leads to the 
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standard bounds on the axion mass, except the possibihty that the axion "starts" at a 
special point seems less viable. We may evaluate the number of axion particles at the initial 
time when the potential is turned on ( defined by the condition m{TQCD) = 'iHiTqcD))) 
na = pa/fTia. Uslug ma{Tf = Vo{T) / ip^ wc may estimate V{Tqcd) ~ T^coi^o/^p^ leading 
to the standard estimate ^ah'^ ~ where h is todays Hubble parameter in units of 100 

km/Mpc/sec. Requiring subcritical fla leads to the standard bound on ipo, ipQ < 
and rria ^ 10~^eV. 

In string theory, natural values of i/jq are approximately lO^^GeV, which, if taken at 
face value, would lead to overclosure of the universe with axions many times over. Two 
possible resolutions have been suggested PD|-P^ to allow our universe to reach its old age 
of today. First, that somehow, perhaps involving some strong coupling string dynamics, 
the low energy effective ipo is some orders of magnitude below lO^^GeV, and the second 
is that some non-standard matter domination epoch, or some late entropy production, has 
occurred in between Tqcd and nucleosynthsis epoch. Our results cannot shed further light 
on this problem, but they do reinforce the need for a resolution. If the resolution of the ipo 
problem requires strongly coupled string theory gi > 1, some of our assumptions should be 
changed but most likely our estimates are still valid, and therefore our results are probably 
qualitatively correct also in that case. Of course, another possible resolution is that the 
model independent axion is not the QCD axion. 

We turn now to the question of energy fluctuations. Since there are large fluctuations in 
the axion fleld (|3.2|) , (p.3|) , we should worry about large energy fluctuations which will cause 
unacceptable deviations from isotropy and homogeneity, affecting either nucleosynthsis or 



the cosmic microwave background. However, as explained in |T8[, these perturbations are 
suppressed. Fluctuations in the axion energy density at a scale i ~ can be computed 
by using the relation 







In previous expressions (■ ■ ■) denotes either vacuum expectation values of operators or sta- 
tistical averages. 

To evaluate (14.31) we need 



(cos (^(0)) cos (^(f))) - (cos (^(0)))^ = e-('^'(°)) [Cosh {{^{0)^{f))) - 1] , (4.4) 
using it we obtain (for the case {ip) = nn, n = 0, ±1, ■ ■ ■) 

J {5pl)^dlnkc:.-V,'e "° , (4.5) 

k 

"'mm 

and, flnally, using (p^) = \Vq we obtain 



1 F ii>i 

i^^^e"'/''^'^, (4.6) 

Pa "00 
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where the upper hmit on the right-hand-side of the previous equation has been changed from 
oo to ks since we take into account only fluctuations produced during the dilaton-driven 
phase. The derivation of eq.(|4.6|) involves some subtleties which we will not discuss. Our 



result agrees with the results of [18 

Ik's 

Because the standard deviation of fluctuations, (yk = J J dlio.k\5-ipk\'^, is much larger 

than the period of the axion potential, ak ^ ipo, energy fluctuations at large wavelength 
are exponentially small, leading to the surprising conclusion that larger field fluctuations 
lead to smaller energy fluctuations. Note that if the spectrum of perturbation is flat as in 
energy perturbations are only power-law suppressed. If, as generally assumed in many 



cases 6tpk < tpQ for all k, then ~ and energy fluctuations actually grow as field 



fluctuations grow. 

Finally, we have not considered any other axion production mechanisms, such as thermal 
production P,^pB|, or the formations of strings, black holes, and other topological objects 
P^ , |25[| which are likely to appear in our model because of the large field fluctuations and 
could result in additional and perhaps dominant axion production leading to a modification 
of our constraints. We hope to discuss these interesting alternatives in the near future. 

The primordial axion spectrum is not thermal, and may consists of a fraction of rela- 
tivistic axions even after their potential is generated. Our understanding of the dynamics 
of the relativistic part of the spectrum is not quite complete because after the axion poten- 
tial is generated the problem becomes an essentially non-linear problem. We believe that a 
better treatment of the relativistic axions is interesting and should be done using numerical 
simulations and tools similar to those used in the theory of topological defects. But we can 
nevertheless reach a few conclusions. First, a necessary condition for a relativistic tail to 
exist after the axion potential is fully developed is that uJsiT ~ ^qcd) > f^a, otherwise 
it can be shown that all modes are non-relativistic. This condition leads to the condition 
zs ^ ^o/Aqcd- If zs > ^o/Aqcd ^ 10^%^o/W^GeV){100MeV/AQCD), then all axions are 
massive. Whether this zs range can be consistent with condition ( |3.7|) depends on i/jq, gi, gs- 
For gi ~ gs the condition becomes IglMp/KqcD ^ zs < iPo/Aqcd, requiring .Ig^Mp > ipo, 
pushing parameters into a relatively narrow region. Modes for which uj/los > AqcdZs/iPo 
are relativistic at T ~ Aqcd- 

The spectrum of axionic perturbations inside the horizon after the generation of the 
potential is quite complicated. A full treatment of these perturbations is outside the scope 
of this paper, and may even result in the conclusion that for the particular case we are 
considering it is not allowed to have any relativistic axions after the potential is generated. 
However, there is a range of frequencies for which we can nevertheless draw definite con- 
clusions. This is the upper end of the frequency range, for which the perturbation is small 
6ijjk < ipoi and relativistic uj > rria. As the universe expands, kinetic energies redshift and 
more axions become non-relativistic. We may evaluate their energy density by calculating 
their number just after the onset of the potential, and, using number conservation, calculate 
their energy density at later times and in particular at matter-radiation equality time. This 

was done in |jl3|, ~ gf (^^^ (^-^^ "L°s'"'^ • Evaluating the axion energy density at 

matter-radiation equality, assuming that all particles have become non-relativistic by then, 
we obtain 
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dhiuj 




2 



3-2^3 




(4.7) 



Since 




10^ 



we see that this region of parameter space gives an uncomfortably large 



energy density, leading to a seemingly favorable region of parameter space zs ^ i^o/ ^qcd- 
However, since our estimates are quite rough, and a small amount of entropy production 
during the evolution of the universe may relax this condition we would not like at this 
moment to completely rule out this interesting possibility. Note that in this case even if 
ipQ ~ IQ^GeV and the axion mass gets pushed towards its upper limit ~ 10~^eV axions 
can provide closure density. This is important for their possible detection. 

In general, for modified spectra, and other relic particles produced by amplification of 
quantum fluctuations during the dilaton-driven phase, it may well be that a fraction of 
relativistic particles remains at teq and therefore it is possible that a single species provides 
simulatneously hot and cold dark matter. 



We have shown that relic axions are produced by amplification of quantum fiuctuations 
with a specific spectrum. In some range of string cosmology model parameters it is predicted 
that most of the energy in our universe today is in the form of cold dark matter axions, 
with suppressed energy density fiuctuations at large wavelengths. Axions could provide 
closure density if their masses lie in the allowed range 10~^eV < rria ^ depending 
on parameters of string cosmology. The spectrum of primordial axions is not thermal, and 
could contain a relativistic tail until quite late times. 
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